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Abstract 

^ . Complex molecules and mesoscopic structures are naturally described by general networks of elementary 

' building blocks and tight-binding is one of the simplest quantum model suitable for studying the physical 
properties arising from the network topology. Despite the simplicity of the model, topological complexity 
' can make the evaluation of the spectrum of the tight-binding Hamiltonian a rather hard task, since the lack 
. of translation invariance rules out such a powerful tool as Fourier transform. In this paper we introduce a 
rigorous analytical technique, based on topological methods, for the exact solution of this problem on branched 
structures. Besides its analytic power, this technique is also a promising engineering tool, helpful in the design 
of netwoks displaying the desired spectral features. 

■ 



I i; buonsaiit@fis.unipr.it 

'"^ ' ^burioniOf is .unipr . it 

^ "^cassiOf is .unipr . it 

o 



1 



1 Introduction: The Role of Topology 



Physical systems like complex molecules and mesoscopic structures consist of a large collection of elementary 
building blocks, atoms or small molecules, connected to give a more or less ordered structure. When a regular 
and periodic pattern is present, as in crystals, it is possible to isolate an elementary cell from which all the 
features of the system stem. At length scales much higher than the cell's size these features depend uniquely on 
the dimensionality of the cell's packing. Periodicity is particularly useful in the analysis of the spectral properties 
of the system: it makes it possible to introduce a (pseudo) wave vector, k, which labels every quantity of interest, 
such as the energy or the Green's functions When, on the other hand, no periodic pattern is recognizable, 
such a simple picture does not apply. In this paper we will develop an analytical method fit for the study of the 
spectral properties of a wide class of discrete structures. We will focus on the simplest quantum model, namely 
the tight-binding model, defined on structures lacking the regular arrangement typical of perfect crystals. In 
order to analyze the features strictly arising from the topology of the arrangement we will simplify as much as 
possible the "internal structure" of the building blocks. The resulting model is simple yet rich, and it displays 
all the qualitative effects due to topology. In this framework we will illustrate an analytical technique, which 
we call bud-reduction, useful for a large class of branched networks. As its name suggests it makes it possible to 
account the spectral features stemming from an entire substructure of a given discrete network by reducing it 
to what we call its bud. This operation is a sort of topological renormalization, in the sense that the presence 
of an entire substructure is accounted by means of an on-site potential of geometrical origin. Moreover it can 
be performed recursively, this meaning that a structure can be reduced to its bud either as a whole, or after 
some of its substructures have already undergone bud-reduction. From a top-down point of view this technique 
allows a geometrical simplification of the original network and a direct understanding of the effect of the network 
topology on its spectral properties. On the other hand it seems particularly well suited for the engineering of 
spectral properties: once the desired spectral structure has been identified it can be obtained by means of an 
appropriate network design. 

The present paper is organized as follows: Section Q briefly deals with the mathematical representation 
of discrete networks, from the topological point of view. A number of useful definitions are also given there. 
Section ^ contains an introduction to the tight-binding model, generalized to the case of an arbitrary network. 
In Section ^ we introduce some mathematical tools which allow the computation of the spectrum of such an 
arbitrary network, independent of its symmetries and regularities. In Section ^ we describe bud-reduction and 
provide numerous examples which clarify the technique and its application. Section ^ consists in our conclusions. 
Section ^ is devoted to some detailed calculations. 

2 Representation of an Arbitrary Discrete Network : Graphs 

When no periodicity is present lattice cannot be used to describe discrete structure, the appropriate choice 
being graphs A graph, Q = {S, A}, consists of a set of points, or sites, S, and a set of links, A, connecting 
points pairwise. In the present paper lowercase italic letters i, j, k,..., will denote the sites of a network, and 
pairs of lowercase letters enclosed in parentheses, {h, k), will denote the link joining the relevant sites. From an 
algebraic point of view a graph is completely described by what is called its adjacency matrix. Every entry of 
this matrix corresponds to a couple of sites of the network, and it equals one if and only if this couple is joined 
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by a link, otherwise it is zero: 

^ _ I 1 if (i, j) is a link of the graph 
^-^ jo otherwise 

Any sequence of consecutive links, {i, h), {h, k), {I, m), (m, j), is referred to as a path on the graph. In the 
following a path starting at site i and ending at site j will be denoted by the symbol Vij. In general such a path 
may reach its ending site, j, also at an intermediate step. When the path never passes through its ending site 
at an intermediate step it will be denoted by the symbol J-ij. The bud-reduction technique we will illustrate in 




Figure 1: Graph of a generic network. Filled circles denote branching sites, i.e. sites which are the unique 
intersection of any two connected subgraphs. Unfilled circles denote non-branching sites. 

Section |5| can be performed at sites of the graph we will call branching sites. A site * of a graph Q is a. branching 
site if it is the only intersection of any two connected subgraphs Bi, B2 oi Q; in mathematical terms 

a = {S, A}, = {S„Ai}, S = SiUS2, A = AiUA2, {*} = SinS2, = Ai n A2 (2) 

Any two subgraphs joined by a branching site are referred to as branches. What makes a branching site peculiar 
is the fact that any path joining sites on different branches is forced to pass through it. Figure ||. illustrates the 
concept of branching site. 
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3 Tight-Binding on a Generic Structure 



The tight-binding approximation [Q] is widely used when deaUng with quantum particles on discrete structures, 
such as solids, molecules, atomic clusters. The Hamiltonian of such a system has the form 

i 

where Vi{r) is the potential energy due to the interaction with the i-th atom of the structure. Assuming that the 
wavefunctions of the Hamiltonian may be expressed as linear combinations of atomic orbitals, that each single 
atom of the structure contributes with a single s-type orbital, and that these orbitals in turn have a non-zero 
overlap only if they come from adjacent atoms, the resulting Hamiltonian has the matrix form 

hJ]^ = tAij + azi5ij (4) 

It consists of an off-diagonal term, called hopping term, and of a diagonal term, called local term. The hopping 
term contains the adjacency matrix describing the topological arrangement of the structure and the so-called 
hopping parameter, t. The local term contains the local parameter and the coordination number, Zi = J2j ^ij^ 
namely the number of nearest neighbours of the i-th site. On a periodic lattice (with homogeneous elementary 
cell) this term is site-independent and it can be dropped without loss of generality, since it simply results in a 
rigid shift of the spectrum. When the geometric structure is non-periodic this cannot be done, and moreover 
there is no natural tool such as the Fourier trasform for the diagonalization of the Hamiltonian. Hamiltonians 
of the form (^) are widely used. They also apply to more general systems than the one described above, in 
which for instance the building blocks at each site are not simple atoms but more complex structures]^, ^. In 
spite of the approximations which it has undergone, the model described above is nevertheless rich and capable 
of displaying the features arising from the topological arrangement of the structure. 



4 Green's Functions and Local Density of States 

A route to the spectrum of the Hamiltonian alternative to direct diagonalization passes through the so-called 
Local Density of States (LDOS), which is related to Green's functions by the formula 

Pi{oj) = --liinlm{Gii{uj + ie)) (5) 

Before going into further detail we point out that most of the quantities we will introduce carry an explicit site 
dependence, which in the case of perfect crystals with simple primitive cell usually drops, due to the spatial 
homogeneity. The Green's function is defined by 

G{uj) = [uj-H]-'^ = [uj-Za-tAr'^ (6) 

where Zij = Zi6ij. By using a complete set of eigenfunctions of the Hamiltonian operator 

H\k)=LOk\k) (7) 

one gets 
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thus 



k 

where \i) is the eigenfunction of the position operator at site i. Recahing that 

hmlm = =p7r (10) 

e-+0 lb ie 

from @ one gets 

p.M=5:j(u;-u;fc)||(i|fc)f (11) 

k 

The last formula sheds light upon the significance of the LDOS at site i: it is a sort of projection of the density 
of states (DOS) on that site. The DOS is recovered by simply summing the LDOS over all sites: 

i k i k 

The normalization of the eigenfunctions of the position operator and the completeness of the Hamiltonian's 
eigenfunctions yield the normalization of LDOS: 

/OO /"OO 
di^ Pi{u:) = Y.\\mf duj5{u^-iVk) = Y.\\^\^)f = ^\i) = ^ (13) 

Simple algebraic manipulations of (^ yield 

1 °° r -1 fc °° 

Gij{uj) = — ^ Uuj - Z a)'^At _ = gjJ2 t'' E fi'i ^i/n 5/^1 ha--- ffhfc_i Afe_i j (14) 

^ ^i"fc=0 k=0 hi,h2,...,/ifc_i 

where 

Qiiu) = (15) 

uj — Zia 

is referred to as atomic Green's function because it is the limit of the Green's function Gij for vanishing hopping 
parameter: 

lim Gi j{uj) = 5i j = gi{uj)5i (16) 

t^o to — Zia 

This can easily be shown by means of (^), simply by observing that uj — H = uj — aZ is a diagonal matrix 
when the hopping term vanishes. Since the chain of adjacency matrix's elements Ai^^ A^-^h^... A}i^_^j in (|l^) 
is non-vanishing only if {z, /ii, j} is a connected /c-step path joining the sites i and j, the Green's 
function Gij{uj) can be viewed as a weighed sum over the paths Vij joining i and j: 

Gijiuj) = gjY^ n Phk (17) 

r^j {h,k)ev^j 

where Phk{t^ o) = tgh{^) is the contribution relevant to the step {hk). Due to the extreme simplicity of the 
present model this quantity is a function only of the site from which the step is taken. In a more refined model 
the hopping parameter would be a link variable, t/i^, and this explains the dependence of phk on both of the 
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sites involved in the step. Note that phk is an adimensional quantity. Thus the LDOS is a local quantity, in the 
sense that it refers to a single site of the structure, but it takes into account the whole structure. This because 
it depends on a Green's function, which is a weighed sum over all possible closed paths passing through the site 
under examination. The function 

r*iM=E n (18) 

appearing in (^) is related to the generating function for a random walk (RW) with probability —agh for the 
step (hk): Tij{u) = Pij{-t/a) 0, |. The function 

7ij{^)=Yi n Phk, (19) 

built using only paths constrained as described in Section ^, is in turn related to the generating function for 
the probability of first arrival: 7ij(u;) = Fij{—t/a). The simple RW is recovered by setting w = 0. The relation 
between Fj j and 7^ j is the same as the one between the generating functions of RW Pij and Fij 

Tij{uj) = 5ij +-fij{uj)Tjj{u;) (20) 

When the last site is the same as the first this reads 

Uu;) = \ (21) 
1 - liM 

where Ti{uj) = Tii{uj) and 7^(0;) = ^ii{uj). 

5 Bud-Reduction 



The analogous of formula (21) in the framework of RW is well known and is due to the factorization properties of 



the generating functions Pa and Fa. Since the proof line for the bud-reduction technique is very similar, in the 



following (21) is proved employing a path technique. First of all we label all the relevant paths in some definite 
though arbitrary sequence for enumerating purposes. Thus let Pj = {'P^"'*}^! be the set of all unrestricted 
paths starting and ending at site i, and let Fj = {Ti^^}n=i be the subset of P consisting of the paths restricted 
as explained in Section ^. Note that both N and N may be possibly infinite. Let txi denote the operation of 
joining two paths. The only requirement for this operation is that the starting site of the second path is the 
same as the ending point of the first one. For instance Vij 1x1 Vjk is a path starting at i and ending at k. Let 

P(.'P)= n (22) 
(hk)ep 

be a function of paths. The following equality holds: 

Pin J M Vjk) = Pi'PiM'Pjk) (23) 

As already pointed out Fj C Pj. Moreover any concatenation of elements in Fj is in Pj. More precisely for 
any V^f G Pj there exists a subset of Fj, {Jf^^f^j^^i, such that pf^^ = j^^^ ex T^f ixi ... JT^^^ Here n is 
nothing but the number of times the path 'P-'^-' reaches site i. Note that the elements of the subset of Fj are 
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not necessarily all distinct. Thus, for any h, pCpI^"^) is the product of of factors of the p{J'ii') kind, each raised 
to some positive integer power n{k). At this point it is easy to understand that the following formal equalities, 
which yield (^), hold: 



N 



k=l 



E 

h=0 



N 
1=1 



(0^ 



E7f 

h=0 



(24) 



Now let us focus on the branching point. Let * be a branching point splitting the graph G into two branches 
as in (^). In the following we are going to show that the evaluation of the Green's function at any site of one of 
the branches can be equivalently performed by pruning out the other branch, provided that the atomic GF at 
the branching point is changed into an appropriate function. It will turn out that the appropriate new atomic 
GF is nothing but the GF of the pruned branch at the branching site. In formula: 



a 



r 



(25) 



where T and B denote respectively the branch containing i, which we call the trunk of the structure, and the 
branch off-springing from the branching site *. As superscripts, these letters denote restriction to the relevant 
substructure. First of all we consider the GF at the branching site *. Since the paths involved in 7^,^, never go 
through * at an intermediate step, they take place entirely either on T or on B. Thus 7** is the sum of two 
terms: 

7*= E Pi'P**)+ E P(^**)=7f + 7f (26) 



From each p{V^^) appearing in ( pq ) it is possible to factor out a term g^. Thus appears in 7^,*, 7^^ and 
only as an overall factor. This means that g^ can be changed into another function, say g.^, through the simple 
multiplication by the factor g^^g^,- For instance 

-1 



l*{^^)\g,=g, = g*i*{i^) (27) 

and similarly for any 7-like function. Such a substitution of the atomic GF can also account for the introduction 
of an on-site potential Vi 
formula (p!5[) would read 



of an on-site potential Vi. In this case the Hamiltonian (^) would have a further diagonal term Vij = 6ij and 



5i(^^) 



ZiQ-Vi 



Now, making use of (^) and (21) it is possible to write 

g* g* 



7f 



1 



7f 



9* 



g* 



1 



7f 



1 -1 



When we restrict the structure to T or B the Green's functions at * read respectively 



r 



Thus, recalling (p] 



becomes 



G 



B 



G 



l-g:^Gtl'^^ 



g* 



7 



r 



G 



r 



g* 



Gl 



1-7 



r 



n -1 



G 



T 



(28) 



(29) 



(30) 
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which is nothing but ( [25| ) in the special case i = *. The proof for the general case can be obtained in the 
same fashion, and it is illustrated in some detail in Section 0. Formula (|25| ) holds for any branching site of the 
trunk. Thus the GFs at the sites of the trunk can be evaluated by pruning off all the branches, provided that 
the atomic GFs at the branching sites are changed into the corresponding GFs of the pruned branches. Since 
they carry all the information relevant to the corresponding branches these new atomic GFs can be viewed as 
their buds. This operation of bud-reduction can be repeated iteratively within the branches and gives rise to 
a topological simplification of the original structure. Of course for an exact result the knowledge of the exact 
Green's functions of the pruned branches is still needed. When these functions are not known a trial expression 
can be plugged in, in order to get approximate results. From a down-top point of view this approach seems 
very well suited to engineering: the effect, in the respect of the LDOS, of appending a secondary structure to a 
known one is simply obtained by plugging in the appropriate atomic GF. Since we will examine only the GFs 
relevant to the same site, from now on we simplify the notation by dropping one of the repeated subscripts. A 
superscript will denote the structure relevant to the GF under examination. The following subsections contain 
a brief survey of possible structures and some examples of bud-reduction. 

5.1 Fundamental Building Blocks 

In this Subsection we apply the technique above illustrated on some simple structures, which can be used as 
building blocks of more complex networks. 



Dimer: a dimer is a simple structure whose graph consists of two sites connected 
by a link. Let gi and gr be the atomic GFs at its two sites. Since there is only 
one restricted path starting and ending at either of its sites, it is easy to show that 
7z = 7r = 9r- Thus, for a dimer 

Gf{gi,gr) = —^ (31) 

Obviously the GF relevant to the other site has the same functional form. It is enough to swap gi and gr- 



Trimer: a trimer is a simple chain whose graph consists of three sites and two links 
connecting two external sites to a central one. Let gi and gc and gr be the atomic 
GFs at its three sites. Since a trimer can be seen as a dimer branch off-springing 
from a dimer trunk we can write: 

Gi{gi,gc,9r) = Gi{gi,Gi{gc,gr)) = - — ,2 ( , \ 

1 - t-'gcigi +gc) 

Gl{9i,gc,9r) = Gfigi,Gfigc,gr)) = —-^^^—^ (32) 

1 - t''9c{9l +9c) 

r<t( \ ndirid.( \ \ 9ri^ — t'^gi9c) 
Gr[9h9c,9r) = Gr[Gr[gi,gc),gr) = --. To— 7 ^ T 

1 - i 9c\9i + 9c) 
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n-mer: iterating the procedure used for the trimer it is possible to deal with a chain 
of length n, for any n. More generally the atomic GF at one end of a n-mer can be 
obtained from the same function relevant to an (/c)-mer, k < n: 

Gt\gi, .■.,gn) = Gf\gi,...,gk-i,G''^~^^^\gk,gk+i,-,gn)) (33) 

In particular, when k = 2 

GS")(<7i,...,5n) = Gf{g,,G^^-'\g2,...,gn-i)) = (34) 

l-VgiG\ '{g2,...,gn-i) 

The GF at an intermediate site can be obtained combining the results for dimers and n-mers as well: 

G^^\gi, ...,gn) = GfiG^igi, gk),G^r'\gk+i, <?n)) VO < < n (35) 



Ring: an n-site ring can be obtained from n-mer by linking the first and the last 
sites. It is equivalent to an infinite linear chain obtained by repeating infinitely 
many times the same sequence {gi, ...,gn} of the ring's atomic GFs. The periodicity 
constraint can be satisfied as follows: the site under examination is identified with 
the central site of a (2n -|- l)-mer, whose external atomic GFs are set to be equal to 
the one we are interested in: 



Gf'^igi, .■.,gn) =9 = Gi^+i^^^(g,52,53, 51,52,53, -^gn^g) (36) 

Note that the periodicity constraint introduces recursion. As it will become clearer below, this recursion forces 
the GF to satisfy a second order algebraic equation. This kind of equation has in general two solutions, of which 
only one is acceptable. The choice criterion is given by relation (16): the GF at every site must coincide with 
the atomic GF at the same site for vanishing hopping parameter. 

Also note that a linear homogeneous chain can be obtained as a ring of length n (n > 3) whose atomic GFs 
are all equal. 



junction: the central site of a trimer can be seen as a 2-fold junction joining the 
two external sites. A 3-fold junction can be simply obtained by plugging a dimer GF 
into the atomic GF of the 2-fold junction site: 



Go (50,51,52,53) = G'*(5i,Gf;(5o,53),52) 



50 



1-^^50 (51 +52 +53) 



(37) 



9 



The GF at one of the three external sites can be obtained from the one at one end of a dimer by plugging the 
GF of the central site of a trimer into the atomic GF of the other end: 



Gl\go,gi,g2,g3) = Gf{gi,Gl{g2,go,g3)) = l-t^ go (52+53] 



51 



1 -i^5o(5i +52 +53) 



(38) 



An n-fold junction, for any n, can be obtained in exactly the same way, simply by iterating the procedure shown 
above. 



Go\go,gi,...,gn) = Gq" ^''^(^'f (5o,5n),5i, •••,5n--i) = ,2 

J- - r 50 Z^j=i 5n 

G;'j(5o,5i,-,5n) = Gf(5i,G("-i)J(5o,52,53,...,5n)=5i(l-t'5iGSj) 



50 



(39) 




Figure 2: Examples of simple bud-reduction. The dark grey bubbles represent the bud atomic GFs. The white 
contoured structures inside them are the relevant branches. Note that the new atomic GF are always evaluated 
at the branching site; a. Trimer. The two sites on the right site are reduced to a bud. The trimer turns into a 
dimer. The GF at the left and right site of this dimer give respectively the GF at the left and at the central site 
of the original trimer; b. 5-junction. The dimer consisting of the central site and site 5 is collapsed to its bud. 
This way it is possible to evaluate the GFs at the origin or at one of the other four peripheral sites by means of 
the relevant GFs for a 4-junction; c. 5-junction. The 4 junction consisting of the origin and of the peripheral 
sites from 1 to 4 is collapsed to its bud. The GF for the dimer allows to evaluate the GF at site 5. 
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5.2 Infinite Chains 



Semi-infinite homogeneous chain (SHC): it is a semi-infinite linear ciiain with 
constant atomic GF, g. Tlie GF at its origin can be evaluated by regarding the 
structure SHC branching off from dimer trunk: 







(40) 



As pointed out while discussing rings, recursion due to periodicity gives rise to a second order equation. The 
requirement limt^o G^^ = g selects the solution 



Gr%) 



1- Vl-(2tg)^ 
2t^g 



25 



1 + ^1 -(2*5)' 



(41) 



The GF at a site different from the origin is evaluated as the same function for one of the ends of an n-mer, 
provided that the atomic GF at that end is changed into the GF at the origin of a SHC. 



Semi-infinite chain (SC): in a semi-infinite chain the origin is a special site; since it 
is the only site with coordination 1, it may have a different atomic GF. Let therefore 
be 5o the atomic GF at the origin and g the atomic GF at any other site. 



(50,5) 



Gf(5o,GS^^(5)) 



1 

90 



1 



{2tgf 



-1 



(42) 



The procedure for a site different from the origin is the same as the one outlined for the SHC. Both in SC and 
SHC cases (^) provides an alternative way: the GF at the k-ih. site is the same as the one for the origin of a 
SC provided that the atomic GF at the origin is the GF for the end of an appropriate /c-mer. 



O — O — O — O — O Homogeneous infinite chain (HIC): let g be the atomic GF at each of the equiva- 
lent sites of an infinite chain. We can obtain the GF at one of these sites by regarding 
it as a branching site where a SHC branch springs off from a SC trunk. Alternatively 
the whole structure may be regarded as a dimer trunk carrying two SHC branches: 

(^^) = Gf^ (gSHC(,),,) = Gf{Gr^{g),Gr^{g)) = ^^^^^^ (43) 

This structure is equivalent to a "ring" composed of two equivalent sites. Thus the same result as in ( ^3|) can 
be obtained from a recursive second order equation for a trimer. The required GF is the GF at the central site 
of a trimer in whose external atomic GFs we plug the required result itself. The central atomic GF is the same 
as in the original structure. 
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Infinite chain with impurity (ICI): if we change the atomic GF of just one site 
of the HIC above into go we obtain an infinite chain with a single impurity. The GF 
at the impurity site can be obtained as follows: 



5.3 Bethe Lattices and Cayley Trees 



1 -1 



(44) 




b. 



c. 



A 



Figure 3: Four fold Bethe lattice (4-Bl); a. The whole structure is shown. The links are drawn with different 
lenghts just for graphical convenience. The shading shows one possible chain substructure which can be chosen 
as the trunk for the bud-reduction; b. Bud-reduced structure: only the trunk is shown, c. Bud-reduced 
structure: the structure of the main buds is shown. They consist of trimers whose external sites' GF are the 
result of a further bud-reduction, d. Bud-reduced structure: the structure of the secondary buds is shown. 
Note that, due to the recursivity of the structure, they are defined in terms of the main buds. The sites at 
which the GFs are evaluated are surrounded by dashed lines. 



An n-fold Bethe lattice (n-Bl) is an infinite network of sites with coordination n, i.e. of n-fold junctions. 
Note that the HIC is nothing but a 2-Bl. Actually a Bl is the simplest looples homogeneous structure after the 
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homogeneous chain. As in a homogeneous chain ah the sites of the n-Bl are equivalent. A possible way to work 
out the GF at one of the n-Bl's sites consists in considering any HIC subgraph of this structure as its trunk. 
Inside each of the n-2 branches off-springing at every site of this trunk one can recognize in turn a SHC trunk, 
dressed with the same atomic GF as the main trunk. Thus the n-Bl can be seen as a HIC, whose atomic GF, J, 
is the one for the central site of a (n-2)-junction; the atomic GF, P for the peripheral sites of this junction are 
in turn the GF for the origin of a SHC, whose atomic GFs once again equal J. Equations ( |45| ) and Figure 
display the bud-reduction described above: 

G^^i = ( = 4-2)J s{g), S{g), ...) S{g) = ( (45) 

" V ' 

n—2 entries 



The combination of the last two equations of (45), together with and (03), give a self-consistency equation 



for J. The GF for the n-Bl is thus worked out by plugging the expression obtained for J into (43), as required 
by the first of (H): 

n — 2 + n^yl — [n — 1)(2 g t)'' 

The details of this calculation are given in Section 0. As stated above a Bethe Lattice is an infinite homogeneous 
structure, and its sites are all equivalent. Its finite version is known as Cayley tree Q. A n-fold Cayley tree of 
radius r ((n, r)-Ct) is a structure invariant under n-fold discrete rotations centered at its central site, which we 
refer to as its origin. Every site whose distance from the origin is lower than the radius has coordination number 
z = n. The peripheral sites, whose distance from the origin equals the radius, have coordination number z = 1. 
A (n, r)-Ct is thus a non-homogeneous structure. Only the sites placed at the same distance from the origin are 
equivalent. Note that, due to the exponential growth, the number of peripheral sites is of the same order as the 
number of internal sites. As illustrated by Fig. ^ the GF at the origin of a (n, r)-Ct can be obtained by means 
of three bud-reduction equations: 

GQ''''^'~^\gn) = Gl\g,Li{gn),Li{gn),...) Li{gn) = G^^'^''\gn,Li+i{gn),Li+i{gn), ■■■) Ly. = gi (47) 

v ' " V ' 

n entries n— 1 entries 

where is the atomic GF at any of the internal sites (z = n), while gi is the atomic GF at any of the peripheral 
sites (z = 1). Note that the second of the (^) sets a recursion relation which is closed by the third equation. 
Also note that imposing the fixed point condition Lj(g„) = Lj_|_i((7„) = L{gn) and plugging the solution into the 
first equation of (^) one recovers the result for the n-fold Bethe lattice, formula (|46|). As illustrated in Figure 
^ , the bud-reduction for the sites other than the origin is a little more complex, due to the lower degree of 
symmetry. 



5.4 Fern Lattices 

A fern is a plant characterized by an evident self-similarity. The secondary structures branching off from the the 
main stalk mimic the whole structure. This resemblance between the branches and the main structure goes on 
up to a certain order /, at which the branches are simple leaves or stalks. An infinite fern lattice has a similar 
structure: at every site of an infinite homogeneous chain b identical secondary structures branch off. These are 
in turn semi-infinite homogeneous chains carrying b identical tertiary structures, and so on. At some order / 
the b substructures branching off from every site are simple semi- infinite homogeneous chains (see Figure ^). 
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Figure 4: Four fold Cayley tree ((4,4)-Ct); a. The whole structure is shown. The links are drawn with different 
lenghts for graphical convenience. The sites with coordination z = 1 are shaded in grey. b. Bud-reduced 
structures for the GF of the structure. From left to right: site at the origin, site at distance i from the origin, 
peripheral site (at distance r from the origin); c. recursion relations for the bud-reduced structures. The sites 
at which the GFs are evaluated are surrounded by dashed lines. 



Thus a fern lattice is characterized by two parameters, the branching number, b, and the order, f. Most of the 
sites of a fern lattice, namely the ones on the branches of highest order, have coordination number z = 2. All 
the other ones have coordination number z = b + 2. Due to translational invariance along the main stalk, the 
sites of a (6, /) fern ((6, /)F) are naturally divided into equivalence classes distinguished by two parameters: the 
order i of the stalk on which they are placed and the distance d from the stalk of order i — 1. Thus i ranges 
from 0, on the main stalk, to /, on the "leaves", whereas — oo <d<ooili = and d > 1 if i > 0. The GF for 
the sites of a {b, f)F can be obtained as 



where gs is the atomic GF at the sites with coordination number z = b + 2, whereas gi is the same quantity for 



{Oi,J{gs,Li)) 



(48) 




b entries b~l entries 



and 




J{9s,Lf) = gi 



(50) 
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Figure 5: Fern lattices. Filled circles: sites with z = 4. Unfilled circles: sites with z = 2; a. (2, 1) fern lattice, 
also known as comb lattice or fishbone; b. (2, 2) fern lattice; some secondary structures are drawn in light grey 
only for graphical convenience. 



It is interesting to note that the fixed point relation Lj+i = Li sets in some sense / = oo. In this case all the 
sites are equivalent and have coordination number z = b + 2. Thus a nBl can be thought of as a (n — 2, oo) fern, 
and indeed the equations (|4^), (|49| ) and (^), together with the fixed point relation, give ( ^6|) once again. 

6 Conclusions 

In this paper we face the general problem of determining the spectral properties of inhomogeneous structures. 
Since the lack of translational invariance rules out a powerful analytical tool such as the Fourier transform, 
and forces a direct-space approach, the spectral properties are recovered through the evaluation of the Green's 
functions, which are explicitly site-dependent quantities. 

An analytical technique, based on the properties of the GFs in presence of branching sites, is developed and 
illustrated by means of a many examples. This technique, which we call bud-reduction, allows to account for 
an entire substructure of the network under examination by means of a sort of on-site potential of geometrical 
origin. It can also be performed telescopically, meaning that the bud-reduced branches may in turn carry 
substructures which have already undergone bud-reduction. This tool yields a topological simplification of the 
structure under examination and allows a better understanding of the effect of each substructure on the spectral 
properties of the system. For this reason it is very promising from the point of view of reverse engineering: the 
insight gained on the effects of different topologies may be a guide in designing networks displaying the desired 
spectral features. 

7 Appendix: Explicit Calculations 

7.1 Bud- Reduction Formula: the General Case. 

In Section |5| we proved formula (|30| ) , which is nothing but bud-reduction formula ( |25| ) in the special situation 
where the site under examination, i, coincides with the branching site, *. We recall that the latter is the unique 



15 



intersection of the two substructure forming the network under examination. We refer to these structures as 
the trunk T and the branch B. The site under examination, i, belongs to T. 

The general case, (|25|), can be easily understood in terms of paths. The function involves paths which 
reach i only after the last step. They can be divided into two sets: the ones which never visit *, which therefore 
take place entirely on T, and the ones which visit * at least once, and can therefore take some steps on B 
(remember that * is the only "door" to B). We denote these paths respectively J^f^ and J^*^. Any path J^*- 
passing through * has the same structure: it is the composition of 

• a path J^i^:f, which takes place entirely on T and never reaches * at an intermediate step. 

• a path "P^^^K which takes place on both T and B and goes through * an unrestricted number of times. 

• a path which takes place only on T and reaches i only after the last step, without going through * 
any more. 



Thus 



E p^^^^-) E E p^^-^^ + E p^\ 



Ep(-^.*)g**5;' Y.pi^*,i) + Y.pi^ti) 



(51) 



T.pi^^,*)T.pi^^,^)\9,=G.. + T.piHi) 



where we used the fact that appears as an overall factor in X^jr, ^ p{^*,i)- If we restrict ourselves to the trunk 
T we similarly get 

' - - - ' (52) 



7, 



r 



J2p(^i,*)J2p(^*^i^\g*=Gj, + J2p(^i 



where G1^, is nothing but the Green's function at * restricted to the trunk. This means that 7^ of the whole 
structure can be obtained through the same function restricted to the trunk, by the substitution G^^^ = 
According to (^0|) this is accomplished by the substitution = Gf^ in GJ^. Thus, since 5* appears only into 

9i 9i 



Gi 



1 - 7i 1 - 7i 



r 



G, 



T 



7»— G?t 



(53) 



which is nothing but (|2 



7.2 Nested Bud-Reduction: Bethe Lattice 

Equations ( |4^ ) explicitly read 



G\ 



nBl 



J 



^l-{2tJY 



J 



1 - (n - 2) t2 ^ 5 



S 



1- Vl-(2tJ)2 
2t2 J 



(54) 
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As we noted in Section 5.3 the last two equations of ([5^) yield a self-consistency equation, 



2{J -g)=g{n-2){l- ^l-{2tJf] 
whose only acceptable solution has the form 



J{9) 



n 



g-{n- 2)^g^l-{n-l){2gty^ 



2 + 2[{n-2)gtf 

The substitution of (^) into the first of ([5^) gives the (HBI) through the use of the following identities: 



(55) 



(56) 



(n{n-2){gtf-Ry = [l + {{n-2)gt) 
n-(n-2)R 



af[n-{n-2)RY 



n-1 



n{n-2){gtf - R (n - 2) + nR 



where R = y/l-{n- l){2gty^. 
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List of Figures 



1 


Graph of a generic network. Filled circles denote branching sites, i.e. sites which are the unique 




intersection of any two connected subgraphs. Unfilled circles denote non-branching sites 


2 


Examples of simple bud-reduction. The dark grey bubbles represent the bud atomic GFs. The 






white contoured structures inside them are the relevant branches. Note that the new atomic GF 






are always evaluated at the branching site; a. Trimer. The two sites on the right site are reduced 






to a bud. The trimer turns into 


a dimer. The GF at the left and right site of this dimer give 






respectively the GF at the left and at the central site of the original trimer; b. 5-j unction. The 






dimer consisting of the central site and site 5 is collapsed to its bud. This way it is possible to 






evaluate the GFs at the origin or 


at one of the other four peripheral sites by means of the relevant 






GFs for a 4-junction; c. 5-junction. The 4 junction consisting of the origin and of the peripheral 






sites from 1 to 4 is collapsed to its bud. The GF for the dimer allows to evaluate the GF at site 5 


3 


Four fold Bethe lattice (4-Bl); a. 


The whole structure is shown. The links are drawn with different 






lenghts just for graphical convenience. The shading shows one possible chain substructure which 






can be chosen as the trunk for the bud-reduction; b. Bud-reduced structure: only the trunk 






is shown, c. Bud-reduced structure: the structure of the main buds is shown. They consist 






of trimers whose external sites' GF are the result of a further bud-reduction, d. Bud-reduced 






structure: the structure of the secondary buds is shown. Note that, due to the recursivity of the 






structure, they are defined in terms of the main buds. The sites at which the GFs are evaluated 






are surrounded by dashed lines. 




4 


Four fold Cayley tree ((4,4)-Ct) 


; a. The whole structure is shown. The links are drawn with 






different lenghts for graphical convenience. The sites with coordination z = 1 are shaded in grey. 






b. Bud-reduced structures for the GF of the structure. From left to right: site at the origin, site 






at distance i from the origin, peripheral site (at distance r from the origin); c. recursion relations 






for the bud-reduced structures. 


The sites at which the GFs are evaluated are surrounded by 






dashed lines 




5 


Fern lattices. Filled circles: sites with z = 4. Unfilled circles: sites with z = 2; a. (2, 1) fern 




lattice, also known as comb lattice or fishbone; b. (2, 2) fern lattice; some secondary structures 






are drawn in light grey only for | 


graphical convenience 
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